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In this paper we investigate the constraints on the cosmoiogicai parameters, especially for the 
equation of state of dynamical dark energy uide, the inflationary parameters n s , a s and r, the total 
neutrino mass ^ m„ and the curvature of universe Qk , using the simulated data of future Planck 
measurement. Firstly we determine the cosmoiogicai parameters with the current observations, 
including ESSENCE (192 sample), WMAP three year (WMAP3), Boomerang-2K2, CBI, VSA, 
ACBAR, SDSS LRG and 2dFGRS, and then take the best-fit model as the fiducial model in our 
following simulations. In the simulations we pay particular attention to the effects of the dynamical 
dark energy in the determination of the cosmoiogicai parameters. Due to this reason, in order to 
make our constraints more robust, we have added the simulated SNAP data in our simulations. Using 
the present data, we find that the Quintom dark energy model is mildly favored, while the ACDM 
model remains a good fit. In the framework of dynamical dark energy model, the constraints on the 
inflationary parameters, ^ m v and Qk become weak, compared with the constraints in the ACDM 
model. Intriguingly, we find that the inflationary models with a "blue" tilt, which are excluded about 
2a in the ACDM model, are well within 2a region with the presence of the dynamics of dark energy. 
The upper limits of neutrino mass are weakened by a factor of 2 (95% C.L.), say, ^ m v < 1.59 eV 
and ^m„ < 1.53 eV for two forms of parametrization of the equation of state of dark energy. 
The flat universe is a good fit to the current data, namely, \Qk\ < 0.03 (95% C.L.). With the 
simulated Planck and SNAP data, the dynamical dark energy model and the ACDM model might 
be distinguished at 4<r confidence level. And the uncertainties of the inflationary parameters, ^ m„ 
and Q.k can be reduced significantly in the framework of the dynamical dark energy model. We also 
constrain the rotation angle Aa, denoting the possible CPT violation, from the simulated Planck 
and CMBpol data and find that our results are much more stringent than the current constraint 
and will be used to verify the CPT symmetry with a higher precision. 

PACS numbers: 98.80.Es; 98.80.Cq 



I. INTRODUCTION 

With the accumulation of observational data, such as the Supernovae Type-la (SN la) [l|, , Cosmic Microwave 
Background (CMB) |, !, !, @ , Large Scale Structure (LSS) Q and so forth, it is possible for us to unveil, despite not 
conclusively for the time being, some enigmas in cosmology, such as the nature of dark energy (DE), the inflation, the 
total neutrino mass ^2 m„, the curvature of our Universe £Ik and even the possible violation of CPT conservation in 
Cosmology [|[. Dark energy, the mysterious source driving the present acceleration of our Universe, has been studied 
widely in the literature since its first discovery in 1998 [9j. However, the nature of DE, encoded in its equation of 
state (EoS) parameter w, remains controversial. Being the simplest candidate of DE and fitting the current data 
well, the Cosmoiogicai Constant (CC), whose EoS remains —1, suffers from the severe theoretical drawbacks such 
as the fine-tunning and coincidence problems 10]. To ameliorate such problems, the dynamical dark energy models 
were proposed. For example, the Quintessence, whose EoS evolves with the cosmic time and satisfies w(z) > — 1 
fill ]. Phantom with w(z) < —1 [l^] and K-essence with w(z) > — 1 or < —1 [l3[. As addressed in literature, recent 
observations mildly favor the DE models with w(z) crossing the cosmoiogicai constant boundary during the evolution 
Ell E3) Il8l [l9L I20I |2~0 ] . Unfortunately, the EoS of the above models cannot realize such "crossing" behavior 
due to the "No-Go" Theorem [2l|, [22|, [23j]. Quintom, whose EoS can cross the cosmoiogicai constant boundary, is 
mildly favored by the observations and has been investigating extensively since its invention jT^, H3| ■ 

Our universe has experienced at least two different stages of accelerated expansion. One is the current acceleration 
driven by dark energy, the other is the inflation in the very early universe [25l [26j . The mechanics of inflation can 
naturally explain the flatness, homogeneity and the isotropy of our Universe. Inflation stretches the primordial density 
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fluctuations and seeds the presently observed large scale structures and cosmic microwave background radiation. In 
2006, the WMAP group claimed that the simple scale-invariant primordial spectrum does not fit well to the Three- 
Year WMAP data [3]. Alternatively, the Harrison-Zel'dovich-Peebles scale invariant (HZ) spectrum (n s = 1, r = 0) 
is disfavored about 2 ~ 3cr. And the large running of the scalar spectral index is still allowed [27[. 

The aforementioned key cosmological questions might be answered by the virtue of the future high precision astro- 
nomical measurements. Especially, the Planck mission of European Space Agency (ESA) will determine the geometry 
and contents of our Universe by measuring the CMB with unprecedented accuracy [28j . Planck will image the full sky 
with sensitivity of AT/T ~ 2 x 10~ 6 , angular resolution to 5' and frequency coverage of 30 — 857 GHz. The angular 
resolution of Planck measurement is three times superior to the current WMAP observation and the noise is lowered 
by an order of magnitude at around 100 GHz. These significant improvements permit much more accurate measure- 
ments of the CMB power spectra, so that Planck has the very power and unique new capabilities to constrain the 
cosmological parameters. In Ref.[28|], the Planck collaboration has done some sensitivity studies of constraining the 
cosmological parameters with simulated Planck data combined with the future SNAP measurement. They investigate 
the dynamics of inflation, neutrino mass, etc. in the framework of the ACDM model and find that with Planck one 
can get much more stringent constraints on the cosmological parameters. 

In our previous works [H, [H, [l9[ we addressed that the determination of the cosmological parameters, such as 
^2 tn>v j and the inflationary parameters, is highly affected by the dynamics of dark energy model due to the 
degeneracies among the EoS of DE and these parameters. Furthermore, dark energy perturbation plays a crucial role 
in the global fit [15j, |29j. Therefore, it is much more fair and reliable to do the error forecasts of the cosmological 
parameters in the framework of dynamical dark energy model rather than assuming a constant w of DE or the ACDM 
model. In this paper, we study the constraints of m w i &k as well as the inflationary parameters in the framework 
of dynamical dark energy models. Using the simulated Planck data, we make a global fit using MCMC method, 
while paying particular attention of the dark energy perturbation in the full parameter space of EoS of dark energy. 
We also stress the role of Planck and CMBpol to detect the possible CPT violation. To obtain the fiducial value of 
parameters for simulation, we firstly constrain these cosmological parameters from the current observations and find 
the best-fit models. 

Our paper is organized as follows: In Section II we describe the method and the current observational datasets we 
used; In Section III we present our method to do the futuristic simulations in detail; Section IV contains our MCMC 
fitting results using the current and future observations and the last section is our conclusion and discussion. 



II. METHOD AND CURRENT OBSERVATIONS 

In our studies, we have modified the publicly available Markov Chain Monte Carlo package CosmoMC 1 [3(| to 
include the dark energy perturbation when the EoS of DE gets across the cosmological constant boundary as we 
illustrate later. We assume the purely adiabatic initial conditions. Our most general parameter space is: 

P = (u>b,oj c , Qfe, 9 S , r, io 0) wi, f v ,n s , log[lQ w A s ],a s ,r) , (1) 

where u>b = Qbh 2 and u> c = il c h 2 are the physical baryon and cold dark matter densities relative to the critical density, 
51^ = I — O m — Ode is the spatial curvature, O s is the ratio (multiplied by 100) of the sound horizon to the angular 
diameter distance at decoupling, r is the optical depth to re-ionization, f„ is the dark matter neutrino fraction at 
present, namely, 

f = = Y,mv (2) 

~ Pdm 93.105 eV n c h 2 ' 1 1 

A s is defined as the amplitude of the primordial spectrum. We parameterize the primordial power spectrum in form 
of 

n s (k) — n s (k s o) + a s In ^— ^ , (3) 

where k s0 is a pivot scale which is arbitrary in principle, here we set A: s o = 0.05Mpc _1 , and a s is a constant 
characterizing the "running" dn s /dlnk of the scalar spectral index, r is the tensor to scalar ratio of the primordial 



1 Available at: |http: / / cosmologist.info / 
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spectrum. The scalar spectral index n s is related to the primordial scalar power spectrum V x (k) by definition: 



, s dVy(k) , l 



Correspondingly, V x (k) is now parameterized as 311 ] : 



k \ a* ( , ( k x x 



ln7> x (fc) = In A s + K(fe s0 ) - 1) In I— J + y (hi ^— j j . (5) 

For dark energy, we choose the commonly used parametrization of the DE equation of state as (32j: 

wde(«) = w a + wi(l - a) , (6) 

where a = 1/(1 + z) is the scale factor and w\ ~ —dw/da characterizes the "running" of the equation of state. In 
left panel of FiglH we divide the (wq, wi) panel into four blocks by lines wq — — 1 and wq + w\ = — 1 as illustrated. 
In the upper right and lower left parts, w(z) is greater or smaller than —1 corresponding to the Quintessence and 
Phantom models respectively. In the other two parts, w(z) can cross the cosmological constant boundary during the 
evolution which can be realized by the Quintom model. The models of Quintom A crosses — 1 from upside down while 
Quintom B crosses from the other direction during the evolution. And the intersecting point denotes the ACDM 
model. However, if one takes the futuristic evolution of EoS into consideration, parts of the region occupied by the 
Quintessence and Phantom will be replaced by Quintom. More explicitly, in the right panel of FigfTJ we redivide the 
parameter space into six parts by the lines wq = — 1, wo + w\ = — 1 and w\ = 0. Part III is for the Quintessence- like 
models, namely, the equation of state remains greater than —1 regardless of the cosmic time, say, w > — 1 for past, 
present and future. Correspondingly, part VI is for the Phantom-like models. Part 1,11, V and IV are all for the 
Quintom- like models. For the models lie within part I and IV, their equation of state has crossed over —1 till now 
while the EoS of the DE models in part II and V will cross —1 in future. 




FIG. 1: Left panel: The parameter space is divide into four parts to distinguish different dark energy models; Right panel: The 
parameter space is divide into six parts including the future behavior of the EoS of dark energy. See text for details. 



Moreover, we also consider another phenomenological parametrization, namely, oscillating Quintom whose EoS 
oscillates with time and is allowed to cross the cosmological constant boundary: 

wde(i) = wo + wi sin(w;2 hi(a)) . (7) 

This oscillating behavior in EoS can lead to the oscillation on the Hubble diagram [3 3| or a recurrent universe which 
can unify the early inflation and the current acceleration of our Universe [13] • In Refs. jl9l. l33l . I3p.l36l some preliminary 
studies have been presented on this kind of dark energy model. From the latest SNIa paper [37J , one can find some 
hints of oscillating behavior in their Fig. 10 where they used a polynomial fitting. Our sine function has the advantage 
of preserving the oscillating feature of the dark energy EoS at high rcdshift measured by the CMB data. For simplicity 
and focusing on the study at lower redshift, we set W2 to be 37r/2 in order to allow the EoS to evolve more than one 
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period within the redshift range of z — to z = 2 where the SNIa data are most robust. We label the above two dark 
energy parameterizations ([5]) and |(7J) as Para I and Para II respectively throughout this paper. 

When using the MCMC global fitting strategy to constrain the cosmological parameters, it is crucial to include 
Dark Energy perturbation. The conservation law of energy reads: 

= , (8) 

where T MI/ is the energy-momentum tensor of dark energy and ";" denotes the covariant differentiation. Working in 
the conformal Newtonian gauge and normally setting the anisotropic stress perturbation of dark energy to be zero, 
one can derive the perturbation equations of dark energy as follows [38] : 

5' = -(1 + w)(9 - 3$') - 3H(c 2 - w)5 - m(w' + 3H(1 + w){c 2 - tu))-r , (9) 

k 

9' = -n(l-3c 2 s )6 + k 2 (^- + ^) . (10) 

However one cannot handle the Dark Energy perturbation when the parameterized EoS crosses —1 based on the 
Quintessence, Phantom, K-essence and other non-crossing dark energy models. By the virtue of Quintom dark energy 
model, the perturbation at the crossing points is continuous, thus we introduce a small positive constant e to divide 
the full range of the allowed value of w into three parts: 1) w > — l + e;2)— 1 + e > w > — 1 — e; and 3) w < — 1 — e. 
Neglecting the entropy perturbation contributions, for the regions 1) and 3) the equation of state does not get across 
— 1 and the perturbations are well defined by solving Eqs. (|9ll0[) . For the case 2), the density perturbation 6 and 
velocity perturbation 9, and the derivatives of 6 and 9 are finite and continuous for the realistic Quintom dark energy 
models. However for the perturbations of the parameterized EoS there is clearly a divergence. In our study for such 
a regime, we match the perturbation in region 2) to the regions 1) and 3) at the boundary and set: 

6' = , 9' = . (11) 

In our numerical calculations we have limited the range to be |e| < 10~ 5 and we find our method is a very good 
approximation to the multi-field Quintom DE model. In our calculation the initial condition we choose is the adiabatic 
perturbations of dark energy, while the isocurvature perturbation of dark energy can be safely neglected [29| . For 
more details of this method we refer the readers to our previous companion papers [H, HI] • 

In our calculations we have taken the total likelihood to be the products of the separate likelihoods (Ci) of CMB, 
LSS and SNIa. In other words defining \\ , = — 21og£j, we get 

Xh,total — Xl,CMB + Xl,LSS + Xl,SNIu • 

If the likelihood function is Gaussian, \\ coincides with the usual definition of x 2 up to an additive constant cor- 
responding to the logarithm of the normalization factor of C. In the computation of CMB we have included the 
three-year WMAP (WMAP3) Temperature- Temperature (TT) and Temperature-Polarization (TE) power spectrum 
with the routine for computing the likelihood supplied by the WMAP team 0, 0, 0, Q as well as the smaller scale 
experiments, including Boomerang-2K2 [HI], CBI 0, VSA [|l[ and ACBAR (43]. For the Large Scale Structure 
information, we have used the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample 0] and 2dFGRS 
[43l ]. To be conservative but more robust, in the fittings to the SDSS LRG sample we have used the first 15 bins 
only, which are supposed to be well within the linear regime. In the calculation of the likelihood from SNIa we have 
marginalized over the nuisance parameter [44]]. The supernova data we use are the recently released ESSENCE (192 
sample) data 1, 2]. Furthermore, we make use of the Hubble Space Telescope (HST) measurement of the Hubble 



parameter Hq = lOOh km s" 1 Mpc -1 45[ by multiplying the likelihood by a Gaussian likelihood function centered 
around h — 0.72 and with a standard deviation a = 0.08. We also impose a weak Gaussian prior on the baryon 
density Q.bh 2 = 0.022 ± 0.002 (1 a) from the Big Bang Nucleosynthesis [4y]. Simultaneously we will also use a cosmic 
age tophat prior as 10 Gyr < to < 20 Gyr. 

For each regular calculation, we run 8 independent chains comprising of 150, 000 — 300, 000 chain elements and 
spend thousands of CPU hours to calculate on a supercomputer. The average acceptance rate is about 40%. We 
test the convergence of the chains by Gelman and Rubin criteria [i?} and find R — 1 is of order 0.01 which is more 
conservative than the recommended value R — 1 < 0.1. 
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III. FUTURE MEASUREMENTS 

When considering the constraints on the cosmological parameters from the simulated data of future CMB 
(PLANCK 2 and CMBpol 3 ) measurements, the fiducial models are obtained by maximizing the likelihood (the best-fit 
model) using the current observations. 

Firstly we derive the likelihood function for a CMB experiment as given by Ref.[48|, |49[. Assuming the CMB 
multipoles are Gaussian distributed, one can obtain the likelihood function as follows: 



exp 



VdetC 

Irn 

where Di m = [af m , af ml af m \ is the data vector of spherical harmonic coefficients which are contributed from the CMB 
signal sim and the experimental noise ni m : af m = s^ + n^, and C is the theoretical data covariance matrix generally 
given by: 

( GJ T Cj E Cj B \ (Cj T + Nj T CJ E Cj B \ 

C = Cj E G EE C EB = Cf E C EE + N EE C EB . (14) 

\ CJ B c EB c BB ) V C I B c f B C ? B + N i BB I 

In this covariance matrix, C^ x denotes the theoretical power spectra and N xx is the noise power spectra which 
can be approximated as: 



N ? X ' = (nfJ n C) = 8xx>el hm A 2 x exp 



f)2 

1(1 + 1) fwh 



In 2 



(15) 



where #f w hm is the full width at half maximum of the Gaussian beam, and Ax is the root mean square of the 
instrumental noise. Non-diagonal noise terms are expected to be zero since the noise contributions from different 
maps are uncorrelated. Due to the global isotropy, the terms Cj B and Cf B are always set to be zero. In our 
calculations we also assume them to be zero except for studying the possible CPT violation in the later section HV CI 
On the other hand, we can estimate the power spectra from the data as follows: 



I Xf Y I 

C XY = \^ \ a lm a lm\ Qg\ 

; ^ 21 + 1 K ' 



So we can obtain the effective X e s'' 



E( 2/ + !)^y h 



Xc 2 ff ^-21n£ = E(2/ + l)/ sky ( ^L +m M + 3j , (17) 



where / s k y denotes the observed fraction of the sky in the real experiments, A is defined as: 

A = Cj T (Cf E C BB - (C EB ) 2 ) + Ci' E (Cf B C EB - Cf E C BB ) + &i B {Qi E C EB - Cf B C EE ) 

+ cJ E (Of B c EB - cJ E c BB ) + c EE (cf T c BB - (Cf B f) + c EB (cJ E CJ B - crc EB ) 

+ cJ B (cJ E c EB - c EE cJ B ) + c EB (cJ E cJ B - crc EB ) + c BB (cf T c EE - {Cf E ) 2 ) , (18) 

and \C\ and |C| denote the determinants of the theoretical and observed data covariance matrices respectively, 

\c\ = crc? E c BB + zcrcrcr - cncrr - cr^r? - cncr? , m 

\C\ = Cj T Cf E C BB + 2Cj E Cj B C EB ~ Cj T (C EB ) 2 - C EE (Cj B ) 2 - C BB (Cj E f . (20) 

The likelihood has been normalized with respect to the maximum likelihood xls = 0, where C XY = C XY . If we set 
the Cj B and Cf B to be zero, the likelihood function will be reduced to the Eq.(17) of Ref.[48(. Furthermore, we can 
obtain the Eq.(9) of Ref.[18| if we ignore the tensor information. 



2 Available at http://sci.esa. int/science-c/www/area/index.cfm?fareaid=17/ 

3 Available at http://universe.gsfc.nasa.gov/program/inflation.html/ 
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TABLE I. Assumed experimental specifications. We use the CMB power spectra only at I < 2500. The noise parameters At 



and A 


p are | 


given in 


units of /iK-arcmin. 




Experiment 


/sky 


^max 


(GHz) 


^fwhm 


A T 


A P 


PLANCK 


0.65 


2500 


100 


9.5' 


6.8 


10.9 








143 


7.1' 


6.0 


11.4 








217 


5.0' 


13.1 


26.7 


CMBpol 


0.65 


2500 


217 


3.0' 


1.0 


1.4 



In some of our simulations we also consider the gravitational lensing effect on the CMB power spectrum. The lensed 
Stokes parameters /, Q and U which specify the intensity and linear polarization of observed CMB are related to 
the unlcnsed Stokes parameters at the last scattering surface (denoted with a tilde) by X(n) = X(n') — X(n + 5n), 
where X denotes /, Q or U and Sn is the angular excursion of the photon as it propagates from the last scattering 
surface until the present. This deflection angle, Sn, is given by the gradient of the lensing potential v<K n )j 

dr r - r " $(rn,r) , (21) 
rr s 

where r is the comoving distance along the line of sight, s denotes the CMB last scattering surface, and \& is the three 
dimensional gravitational potential (50 . l5~fl ] . 

The important feature is that the gravitational lensing can mix E and B modes (Hoj . If we assume that there is 
only unlensed E type polarization and the unlensed C BB = in the last scattering surface, the gravitational lensing 
will generate B type polarization in the observed field, Cf B ^ 0. The information from the gravitational lensing is 
added through the power spectrum for the lensing potential Cf^ and the correlation to the temperature Cj^: 

(4l4m>) = (Cf 4, + Nff)S w 6 mm> , {a24 m ,) = (Cf * + N^)5 w S mm , , (22) 

which can be computed numerically in the linear theory using CAMB 4 [52| . In our analysis we use the unlensed power 
spectra, Cf T , Cf E , Cf E , and , Cf . We do not use the lensed power spectra to avoid the complication of the 
correlation in their errors between different I values and with the error in Cj and Cf 4 [H, H3] . For errors on Cf^ we 
follow the Ref. 55]. We use the publicly available code 5 [49| to simulate the mock CMB power spectra of our fiducial 
models. In Table I we list the assumed experimental specifications of the future Planck and CMBpol measurements 
and neglect the foreground contamination. 

To make our constraints more robust, we add the simulated SNAP data to do all the simulations throughout this 
paper. The projected satellite SNAP 6 (Supernova/Acceleration Probe) would be a space-based telescope with a one 
square degree field of view with 10 9 pixels. It aims to increase the discovery rate for SNIa to about 2000 per year. 
The simulated SNIa data distribution is taken from Refs.JH, [13] • As for the error, we follow the Ref. (5(| which takes 
the magnitude dispersion 0.15 and the systematic error cr sys = 0.02 x z/1.7, and the whole error for each data is as 
follows: 



/ 15 2 

O-magOi) = \ a\ Uzi) + — , (23) 



where m is the number of supernova in the i'th redshift bin. 



IV. RESULTS 



In this section we show our global fitting results of the cosmological parameters and focus on the dark energy 
parameters, inflationary parameters, space-time curvature, total neutrino mass and the rotation angle denoting the 
possible CPT violation respectively. 



4 Available at |http:/ /camb.info/| 

5 Available at |http://lappweb.in2p3.fr/~perotto/FUTURCMB/home.html/| 

6 SNAP is one of the several candidates emission concepts for the Joint Dark Energy Mission (JDEM). Available at http:/ /snap. lbl.gov 
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TABLE II. Constraints on the EoS of dark energy and some background parameters from the current observations and the 
future simulations. Note that Para I and Para II represent wue(ci) = u>o + Wi(l — a) and uide(i) = wo + Wi sin(37r/2 ln(o)) 
respectively. For the current constraints we have shown the mean values l,2cr (Mean) and the best fit results together. And 
we also list the standard deviation (SD) of these parameters based on the future simulations. 





ACDM 


Para I 


Para II 




Current 


Future 


Current 


Future 


Current 




Best Fit Mean 


SD 


Best Fit Mean 


SD 


Best Fit Mean 


Wo 


-1 -1 


-1.16 —1.03+S il + oll °- 045 


-0.898 -0.9811°;™'.^ 


Wl 


0- 


0.968 0.405i° ;^I°;™ 0.11 


0.047 -0.068±° a j^r° 2 S 


f«DE 


0.760 0.762tomlto°oll °- 0043 


0.756 0.760±ooig±oo35 0.00 64 


0.765 0.76^ZltlZl 


Ho 


73.1 73.3±i;?±l;| °- 44 


70.3 71.2I^t£| 0.76 


72.0 72.2 + ^ +s 6 ;° 



A. Equation of State of Dark Energy 

To study the dynamics of dark energy, we parameterize our universe as follows: 

P= (Lj b ,u c ,Q s ,T,w , Wl ,n s ,log[W 10 A s ]) . (24) 

Our main results of dark energy parameters are summarized in Table II. Besides the two parameterizations Para I 
and Para II, we also investigate the ACDM model for comparison. In addition, we present the future constraints for 
the ACDM model and Para I using the simulated Planck and SNAP data as introduced above. Marginalized over 
other cosmological parameters, in Table II we list the constraints on the dark energy parameters as well as the Hubble 
constant in different dark energy models. 




-1.3 -1.1 -0.9 -0.7 -0.5 -2 -1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 



FIG. 2: Constraints on the dark energy parameters wo and iui from the combination of current observations (Black Solid Lines) 
and the future simulation data (Red Dashed Lines) respectively. The left panel is for Para I: wue(cl) = »o + Wi (1 — a). And 
the right panel is for Para II: wde(«) = wo + w i sin(^ log(a)). The two blue dotted lines in the (too, wi) panel distinguish the 
dark energy models and their intersecting point denotes the ACDM model. 

In Fig|5]we illustrate the constraints on the dark energy parameters wq and w\ of two parameterizations. From the 
current observations we find that wq = —1.03 ± 0.15, w\ = 0.405+0 587 for Para I and the Quintom scenario, where 
w(z) can cross the cosmological constant boundary during the evolution, is mildly favored. Using the current data, 
we find that the best fit model is located in the Quintom A region while the ACDM, denoted by the intersect of 
two straight lines, lies at the edge of la contour. The one dimensional constraint on the evolution of w(a) from the 
current data is shown in FigJ3] This behavior can be found more obviously from the best fit model. However, current 
data can not distinguish different dark energy models decisively, namely, the variance of wq and wi are too large to 
distinguish dynamical dark energy models from the ACDM model. The ACDM model is still a good fit right now. 

In order to distinguish different Dark Energy models we consider the future measurements Planck and SNAP. The 
fiducial model we choose is the best fit model from the current constraints of Para I. We show the 68% and 95% 
confidence level contours (Red Dashed lines) on the left panel of FigJ^l As we expected, the constraints from the 
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0.0 0.2 0.4 0.6 0.S 1.0 12 1.4 1.6 l.S 2.0 

Z 



FIG. 3: Constraints on woe(<i) = Wq + — o) from the current observations. Median value (central red solid line), 68% 
(inner, dark shaded area) and 95% (outer, light shaded area) intervals are shown. The blue dashed lines denote the cosmological 
constant boundary. 

TABLE III. Constraints on some cosmological parameters n 3 , a s , r, fit, and ^ m v from the current observations and the 
future simulations. We have shown the mean 1,2a errors (Mean) for the current constraints and the standard deviation (SD) 
of these parameters based on the future simulations. For the weakly constrained parameters we quote the 95% upper limit 

instead. 





ACDM 


Para I 


Para II 




Mean SD 


Mean SD 


Mean 


n s 


053+0. 014+0. 028 Q nno 
U.9OO_ .013-0 .026 U.UUO 


0.965±° ;°oi£S 0.0037 


n QCO+0. 016+0. 036 
u - yDZ -0. 017-0. 031 


100 x a s 


o 7 r+2.19+4.24 n ro 
— 2.21 — 4.23 U.OO 


q qo+2.52+4.80 n PC 
o -°°-2.50-4.76 U - JJ 


q qp+2.37+4.86 
0.»0_ 2 39 _4 72 


r 


< 0.231 (95%) < 0.055 (95%) 


< 0.392 (95%) < 0.074 (95%) 


< 0.356 (95%) 


too x n k 


-0.873±^™ 0.289 


-0.20lti; 2 6 9 ± 2 2 .S 1.05 


-0.593l™: 5 5 7 




< 0.958 (95%) 0.077 


< 1.59 (95%) 0.179 


< 1.53 (95%) 



simulated data are much more stringent than the current constraints. By the virtue of Planck and SNAP data, we 
see that the standard deviations of wo and w\ are a = 0.045 and a = 0.11 respectively, which are reduced by a factor 
of 3.33 and 5.2. The Quintom model and the ACDM model might be distinguished at around 4er confidence level. 

For Para II, the mean values from the current observations are u>o = — 0.981^q'34q, wi — — O.OGS^QggJ which still 
support the Quintom scenario despite of the weak significance. In right panel of Fig|2l we see that the Quintom 
models occupies the most of the contour while the ACDM model still lies well within la contour. 

B. Other Cosmological Parameters 

The dynamics of dark energy can have profound effects on the determination of other cosmological parameters, such 
as the inflationary parameters (n s , a s , r), the total neutrino mass ^2m u as well as the curvature of space-time fi^, 
due to the well-known degeneracies among these parameters. In this subsection, we measure the above parameters in 
the framework of dynamical dark energy models. 

1. Inflationary Models 

The current acceleration and the inflation, the two stages of accelerated expansion of our universe, might have 
some deep relationship albeit the significant difference of energy scale between them. Some efforts have been made to 
unify these two expansion epoches, such as the quintessential inflation (58|. Moreover, the isocurvature perturbations 
in dark energy sector generated during inflation may give rise to the suppression of power of CMB at large scale, 
which can be mimicked by suppressed primordial spectrum [5t| . This means different dynamics of the dark energy 
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and inflation can lead to similar effects on observations and studying such degeneracies might unveil the possible 
connections between dark energy and inflation. 




0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 4: ID current constraints on the parameters n s , a s and r based on the different dark energy models: ACDM (black solid 
line), wde = Wo + wi(l — a) (red dashed line) and wde = wo + w i sin(^ log(a)) (blue dotted line). 

In this section, we constrain the inflationary parameters in the framework of dynamical dark energy using the 
current and simulated datasets. We sample in the following 10 dimensional parameter space using MCMC algorithm: 

P = (u} b ,uj c ,e s ,T,w ,w 1 ,n s ,\og[10 w A s ],a s ,r) . (25) 

It's noteworthy that we do not constrain a s and r simultaneously in our global fittings. From Table III, we can see 
the effects of dynamical dark energy on the determination of the inflationary parameters. Again, we give the fitting 
results for Para I, Para II and the ACDM model for comparison. We find that the constraints for the spectral index 
n s , the running a s and the tensor-to-scalar ratio r have been weaken with the presence of dynamics of dark energy. 
Quantitatively, the 2a constraints of n s , a s and r can be relaxed by roughly 36%, 13% and 70% respectively. This 
can be seen from the one dimensional distribution plot of FigJU 

The WMAP group found that the scale invariant primordial spectrum and the inflation models with n s > 1 is 
disfavored at almost the 3a level. Our result is in good agreement with them, n s — 0.953^Q g}g, based on the ACDM 
model. However, we find that the mean value of n s moves toward to the "blue" spectral in the framework of dynamical 
dark energy model, n s — 0.965 ± 0.017. From the future data we find that the standard deviation of n s can be shrink 
to be a — 0.003 and the scale invariant spectrum will be verified with much higher confidence level. 

In the framework of dynamical dark energy model, from Figdl we find that the 95% upper limit of r can be relax 
from r < 0.231 to r < 0.392. The degeneracy may be due to the reason that the tensor fluctuation and the dark 
energy component mostly affect the large scale (low multipoles) power spectrum of CMB. In the two dimensional plot 
of FiglSl we find that the Harrison-ZePdovich-Peebles scale invariant (HZ) spectrum (n s = 1, r = 0) is disfavored 
about 2 ~ 3a in ACDM model. However, this spectrum can be allowed with the presence of the dynamics of dark 
energy. The single slow-rolling scalar field with potential V(4>) ~ m 2 2 , which predicts (n s ,r) = (1 — 2/JV, 8/JV), 
is well within la region, while another single slow-rolling scalar field with potential V(<fi) ~ A</> 4 , which predicts 
(n s ,r) = (1 — 3/iV, 12/7V), is excluded about 2a in the ACDM model. Interestingly many hybrid inflation models, 
excluded in the ACDM model, revive in the framework of dynamical dark energy model as illustrated in FigfSl 

Another feature of WMAP data, both for WMAP1 [60| and WMAP3 @, H3, is the large running of the scalar 
primordia spectrum index a s . Our result shows that the large running is favored more than la, a s — —0.038 ± 0.022. 
In FigU we find that the dynamical dark energy models enlarge the error of a s slightly and do not affect the mean 
value obviously. 
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FIG. 5: 68% and 95% constraints on the panel (n s , r) based on the different Dark Energy models: ACDM (black solid line), 
hide = wo + Wi(l — o) (red dashed line) and woe = wo + wi sin(^ log(a)) (blue dotted line). The two solid green lines delimit 
the three classes of inflation models, namely, small-field, large-field and hybrid models. The blue points are predicted by m 2 </> 2 
model and \<t> 4 model respectively. These predictions assume that the number of e-foldings, N, is 50 — 60 for m 2 (f> 2 model and 
64 for X(f> 4 model. The magenta dash-dotted lines denote the 1, 2a contours obtained from the future simulated data. 



Given the large uncertainties in the constraints of inflationary parameters from the current observations, different 
inflation models cannot be distinguished conclusively. Yet, the constraints from the future Planck measurement can 
make this possible. From our simulation results in Table III, we find that the error bars of inflationary parameters 
can be reduced by about a factor of 5. This dramatic improvement will play a crucial role in the study of dynamics of 
inflation and can also shed light on the investigate of the dynamical dark energy model due to the correlations among 
inflationary and dark energy parameters. 



2. Curvature of Universe 



Dark energy and the curvature, Qk — 1 — ~ ^de, are dominant factors in determining the fate of our Universe. 
Further, DE parameters and £Ik are correlated. This is expected since £Ik and dark energy can contribute to the 
luminosity distance via: 



E(z) 



d L (z) 
H 



l + z 



zsmn 



dz' 
o W) 



^On(l + Z) 3 + r>DECX P ^3jf 1 X^ dz' 



+ n K (i + z) 2 , 



(26) 
(27) 



where sinn(y \k\x)/w \k\ — sin(x), x, sinh(x) if k < 0, k = and k > 0. In addition, can modify the angular 
diameter distance to last scattering surface, which leaves imprints on the CMB power spectrum. 

We concentrate on the determination of Qk m dynamical dark energy models using current and simulated data. 
Our parameter space is: 



P EE (LJ b ,UJ c ,Qk,Q S ,T,W0,Wl,n s ,log[10 A s }) 



(28) 



From Table III and Figl6l we see our universe is very close to flatness, namely, the absolute value of space-time 
curvature \Clk\ is smaller than 0.025 in the ACDM model, 0.028 for Para I and 0.032 for Para II. The dynamics of 
dark energy weakens the constraint of \Qk\ due to the well-known correlation among f2^ and dark energy parameters. 
This correlation plays a crucial role in the reconstruction of equation of state of dark energy [61| . By the simulated 
data, we are able to detect the curvature more accurately. 
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FIG. 6: ID current constraints on the parameters fl*, based on the different dark energy models: ACDM (black solid line), 
hide = too +m>i(1 — a) (red dashed line) and hide = wo + w\ sin(^ log(a)) (blue dotted line). 

3. Neutrino Mass 

Detecting the absolute mass of neutrino is another challenge of modern physics. The cosmological observations can 
obtain upper limits of the absolute neutrino mass. For background evolution, neutrino masses, albeit small, contribute 
to the cosmic energy budget and modify the epoch of matter-radiation equality, angular diameter distance to the last 
scattering surface and other related physical quantities. For the evolution of perturbation, neutrino becomes non- 
relativistic at late time thus they damp the perturbation within their free streaming scale. Thus the matter power 
spectrum can be suppressed by roughly AP/P ~ — 8£l v /Q m [62]. As a result, neutrino can leave imprints on the 
cosmological observations, such as CMB and matter power spectrum. On the other hand, the evolution of dark energy 
can also affect the evolution of background and perturbation, which mimics the behavior of neutrino to some extent. 
This leads to an obvious degeneracy among dark energy parameters and the neutrino mass. 

The degeneracy between dark energy with constant equation of state and neutrino mass has been studied in the 
literature [1, [63[ . In this section, we update our previous results to study the upper limits of neutrino mass with the 
presence of dynamical dark energy [18j and investigate the degeneracy between dynamical dark energy and neutrino 
mass with current cosmological observations as well as with the future simulated data. 
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J 
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FIG. 7: ID current constraints on the parameters ^ m„ based on the different dark energy models: ACDM (black solid line), 
woe = Wo + wi(l — a) (red dashed line) and hide = wo + mi sin(4^ log(a)) (blue dotted line). 

We concentrate on the determination of ^2 m v in dynamical dark energy models using current and simulated data. 
Our parameter space is: 

P = (wt, lo Ci e s , t, w , wi,f v , n„ logllO 10 ^]) . (29) 

In the last row of Table III, one can read 95% C.L. neutrino mass limits derived from the current observations as well 
as the simulated data of Planck and SNAP in the ACDM and dynamical dark energy models. In the ACDM model, 
the limit of neutrino mass we get, ^2 m v < 0.958 eV (95%), is consistent with Tegmark's result 0- For the dynamical 
dark energy model, the limit can be relaxed to m v < 1-59 eV (95%) and m v < 1-53 eV (95%) obviously, due to 
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the deg eneracy between the dark energy parameters and the neutrino mass from the geometric feature of our Universe 
[II If. In FigEI we illustrate this effect with current astronomical data. 

With the simulated data, we obtain the two tail posterior distribution due to the nonzero fiducial value of neutrino 
mass. The standard deviation will be greatly tightened to be 0.077 eV. With the presence of dynamical dark energy, 
the standard deviation can be relaxed by a factor of 2.3 using the simulated data. We might distinguish the normal 
hierarchy (^nit, ~ 0.05 eV) from the inverted hierarchy (^2,m u ~ 0.1 eV) using the future Planck measurement. 



C. Cosmological CPT Violation 



The CPT symmetry which has been proved to be exact within the framework of standard model of particle physics 
and Einstein gravity could be violated dynamically during the evolution of the universe The detection of CPT 
violation will reveal new physics beyond the standard model. In our previous work we studied the cosmolo gica l CPT 
violation in the photon sector. We introduce a Chcrn-Simons term in the effective Lagrangian of the form [651 ] : 

AC = ~P„A V F»» , (30) 

where p^ is a four- vector and F^ v = (l/2)e A " yp<T F pcr is the dual of the electromagnetic tensor. This action is gauge 
invariant if p^ is a constant and homogeneous vector or the gradient of a scalar. It violates Lorentz and CPT 
symmetries if the background value of p^ is nonzero. In the scenario of quintessential baryo-/leptogenesis [66L l67| the 
four-vector p^ is in the form of the derivative of the quintessence scalar, d^4>. During the evolution of quintessence, 
the time com pon ent of d^cj) does not vanish, which causes CPT violation. In the scenario of gravitational baryo- 
/leptogenesis [681 , |69| , p^ is the gradient of a function of Ricci scalar R. 




FIG. 8: ID current constraints on the rotation angle Aa from the simulated Planck (Black Solid line) and CMBpol (Red 
Dahsed line) data. 



The interaction in Eq. ([30|) violates also the P and CP symmetries, as long as po does not vanish [70j| . It leads to a 
rotation of the polarization vector of electromagnetic waves when they are propagating over cosmological distances. 
This effect is known as "cosmological birefringence" . The polarization vector of each photon is rotated by an angle 
Aa and one would modify the power spectra of TE, EE, BB, TB and EB at the last scattering surface as: 

C{ TB = Cr B sin(2Aa), (31) 

C[ EB = \{Cr ~ Cr>H^) , (32) 

C[ TE = Cj E cos(2Aa) , (33) 
C[ EE = Cf B cos 2 (2Aa) + Cf s sin 2 (2Aa) , (34) 
C[ BB = Cf B cos 2 (2Aa)+Cf B sin 2 (2Aa) , (35) 

where the primed quantities are rotated. In Ref.Q we have performed a global fit and found that a nonzero rotation 
angle of the photons is mildly favored, Aa = —6.0 ± 4.0 deg, using the WMAP3 (without the information of TB and 
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EB power spectra) and full data of Boomerang-2K2 7 . 

In FigJSl by the simulated data with Planck and CMBpol accuracy, we find that the standard deviation of the 
rotation angle will be greatly tightened to a = 0.057 deg for PLANCK and a = 2.57 x 10~ 3 deg for CMBpol. These 
results are much more stringent than the current constraint and will be used to verify the CPT symmetry with a 
higher precision [73j]. 

V. SUMMARY 

Since the mystery of our Universe encodes in the cosmological parameters, constraining these parameters with the 
latest observational data and doing error forecasts with the simulated futuristic data can lead us to understand the 
Nature. In this paper, we focus on the dynamics of dark energy in light of current and simulated Planck and SNAP 
data and then constrain the inflationary parameters, total neutrino mass and curvature of space-time in the framework 
of dynamical dark energy model. In addition, we investigate the rotation angle Aa, the possible signature of CPT 
violation, by simulated Planck and CMBpol data. 

Parameterizing the EoS of dark energy in two forms as Eqs. (|6l7p . we find that the Quintom model, whose EoS 
crosses —1 during the evolution, is mildly favored by latest observations, albeit the ACDM model remains a good 
fit. Using the simulated Planck data complimented by SNAP data, we find that the variance of the dark energy 
parameters in Eq.Q decreases dramatically, namely, a(wo) and a(wi) can be reduced by a factor of 3.33 and 5.2 
respectively. Given the current central value, this means that Planck can distinguish dynamical dark energy from the 
ACDM model at around 4a confidence level. 

Since the dynamics of dark energy greatly affects the determination of other cosmological parameters, we constrain 
the inflationary parameters, total neutrino mass ^ m v and curvature of space-time Q.k with the presence of dynamics 
of dark energy using current and simulated observational data. We find that the dynamics of dark energy generally 
weakens the determination of other cosmological parameters. For instance, we find that the inflation models with a 
"blue" tilt (n s > 1), which is strongly disfavored in the ACDM model, are now well within 2a region in the framework 
of dynamical dark energy model. This is something intriguing in that Hybrid Inflation Models have been revived due 
to the dynamical dark energy, due to the obviously enlarged parameter space of (n s ,r) in the framework of dynamical 
dark energy model. These discoveries can lead us to further study the dynamics of inflation, dark energy and the 
relationship between them. With Planck, the uncertainties of inflationary parameters, including n s , a s , r, can be 
roughly reduced by a factor of 5. This significant improvement make it possible to distinguish the inflationary models 
decidedly. 

The presence of dark energy can relax the upper limit of neutrino mass by a factor of 2. This sheds new light 
on the research of neutrino physics, such as the scenario of mass varying neutrino [zl, [75[. By Planck, we can rise 
the precision of measurement of this mass limit by a factor of 12, namely, the standard deviation of neutrino mass 
can be shrunk to 0.077 eV for the ACDM model and 0.179 eV for the dynamical dark energy model (Para I). This 
measurement might make us understand the Nature of physics conclusively. Current data imply that our Universe is 
very close to flatness, say, \Qk\ < 0.03 (95% C.L.). By Planck, we can reduce a (fix) by a factor of 2.67(ACDM). 
Such improvement helps us to reveal the Nature of space-time. 

The global symmetry of CPT plays a critical role in understanding the fundamental physics. In our previous works 
we have found some hints of CPT violation encoded in the rotation angel of polarization vector of photons Aa. With 
Planck and CMBpol, we can test the CPT symmetry at a unprecedent accuracy. 
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7 The WMAP group did not release the TB and EB data when we prepared our paper Q at that time. We had to set the TB and EB 
of power spectra of WMAP3 data to be zero, C ; TS = and Cf B = 0. Recently Cabella et al. 17111 performed a wavelet analysis of 
the temperature and polarization maps of the CMB delivered by the WMAP experiment which includes the information of TB and 
EB power spectra. They obtained a limit on the CMB photon rotation angle Aa = —2.5 ± 3.0 deg. Right now the WMAP group has 
released their results of TB and EB power spectra. We plan to combine the full data of WMAP3 (including the information of TB and 
EB power spectra) and BOOMERANG-2K2 to reanalyze the rotation angle in the future [72| . 
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